Analytical solutions for two-level systems with damping 
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A method is proposed to transform any analytic solution of the Bloch equation into an analytic 
solution of the Landau-Lifshitz-Gilbert equation. This allows for the analytical description of the 
dynamics of a two level system with damping. This method shows that damping turns the linear 
Schrodinger equation of a two-level system into a nonlinear Schrodinger equation. As applications, 
it is shown that damping has a relatively mild influence on self-induced transparency but destroys 
dynamical localization. 



Two-level systems have become almost ubiquitous in 
modern physics. They are found for instance in laser 
physics, magnetic resonance spectroscopies, quantum 
computers, quantum teleportation and optoelectronics. 
The state of a two-level system can be described by 
an effective moment M and its dynamics by the equa- 
tion M = — 7M X B, where 7 is the gyromagnetic fac- 
tor and B is an external time-dependent field. We fol- 
low the standard (and inappropriate) custom of calling 
M = -7M X B the "Bloch equation" . 

The main drawback of the Bloch equation is its absence 
of damping. There are many phenomenological models 
of damping for two-level systems. When decoherence is 
small, the system remains in a pure state, the length of 
M is constant and damping is taken into account by the 
so-called Landau-Lifshitz-Gilbert (LLG) equationP,0,[l| 
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where M = |M|, M = dM/dt and a > 0. It was shown 
that the LLG equation provides a realistic model of ferro- 
magnetic resonance, micromagnetics, spin- valve dynam- 
ics [J|, the magnetism of thin films Q and nanomagnets 
d|, the dynamics of domain walls in various geometries 
3. Cheng and coll. recently proved that the LLG equa- 
tions could be derived from a microscopic model [Sj] . 

The linear nature of the Bloch equation allowed for 
the discovery of many analytic solutions. For example, a 
recent work identified 26 families of solutions that can be 
expressed in terms of special functions [9| and significant 
steps toward the general solution were carried out p^ . 
Even when analytic solutions are not available, powerful 
analytic approximation methods exist [111, [H, • 

By contrast, very few solutions of the physically more 
accurate LLG equation are known. In this paper, we de- 
scribe a method by which any analytic solution of the 
Bloch equation can be transformed into an analytic solu- 
tion of the LLG equation. Similarly, any analytic approx- 
imate solution of the Bloch equation is transformed into 
an analytic approximate solution of the LLG equation. 
Then, we show that this transformation turns the linear 
Schrodinger equation for a two-level system into a non 
linear Schrodinger equation. Finally, we investigate the 
influence of damping on self-induced transparency and 
dynamical localization. 



We describe now the transformation from the solu- 
tion of the Bloch equation to the solution of the LLG 
equation. Consider a solution of the Bloch equation 
M(7) = — 7M(7) X B, where the dependence of AI on the 
gyromagnetic factor was written explicitly and where B 
is a real function of t. Assume now that M(7) is an an- 
alytic function of 7. This allows us to define N = M(7), 
where 7 — 7/(1 — ioi). This gives us 



N = M(7) = -7M(7) X B 
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-N X B. 



The equation of motion implies that 
not depend on t. We define now 



^ = 



A^x + iNy 



E.Mf does 



(2) 



If we calculate the derivative of ^ with respect to t, taking 
account of the fact that M does not depend on time, we 
find 



+ iNy (A^:, -I- iNy)N, 



(M + A^, 



(3) 



If we substitute the equation of motion for N, we can 
check that ^ satisfies 



2{1 ~ia) 



(B-e' + 2B,e-B+), (4) 



where = ± iBy. The complex function ^ is used 
to define a real vector M' by 



ML 



so that we still have 



M, 



M, 



M + M'' 



(5) 
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but Mj' are now real. The question is: what is the equa- 
tion satisfied by M'? Let us calculate the derivative of 
M' with respect to t. We have 
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M, 



If we express ^ and ^* through equation (g]) and its con- 
jugate, we obtain M' in terms of ^ and If we replace 
them by equation ([5]) and its conjugate, we obtain, after 
a lengthy but straightforward calculation. 
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x{ByM'^M' - B,M'^ + B.M'^M', - S.A/f ) 
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x(i?,M;M^ - i?,Mf + B^M'yM'^ - ByM'^^) 



ML = -- 



x{B^M',M'^ - B.M'J" + ByM'.M' - B,M'^) 



This can be rewritten 
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rM' X B 



70; 



(1 + a2)M 



M' X (M' X B) 



We recognize the Landau-Lifshitz equation, in a form 
equivalent to the LLG equation. To show the equiv- 
alence, use the LLG equation to derive M' x M' = 
— 7M' x (M' X B) — aA/M', introduce this expression 
in the LLG equation and solve for M'. Thus, we have 
transformed a solution M of the Bloch equation into a 
solution M' of the LLG equation ([T]) . 

We discuss now the simplest example of this transfor- 
mation, to show how the method works in practice. We 
consider a constant external magnetic field along the z 
axis, i.e. B(t) — (0,0, i?o)- The solution of the Bloch 
equation is: 

M^ = M sin 6*0 cos{nt + 0o), 
My = Msin6'osin(fit-|- (/)o). 



M cos 6^0, 



where Af, and (f>Q are constants and Vi = jBq. M is 
obviously an analytic function of 7 and we can substitute 
7/(1 — ia) for 7 in M to define the vector N. This gives 



with n' = + a^) and we recover the well-known 

solution of the LLG equation 
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sin^o cos(il'i + 4>o) 
cosh{afl't) + cos 6*0 sinh{ail.'t) ' 

sin^o sin(f7'i -I- (j)o) 
cosh(ail't) + cos 6*0 sinh(an'<:) ' 
cos ^0 cosh(ari't) + smh{a^l't) 
cosh{afl't) + cos 6*0 sinh(ari't) 



If > 0, the equilibrium magnetization is M' = 
(0,0, Af); if f7 < 0, it is M' = (0,0,-Af). As expected, 
damping transforms a precession dynamics into a motion 
towards an equilibrium state. 

It is convenient to determine directly the influence of 
damping on the dynamics of the two-level system de- 
scribed in the Schrodinger or Heisenberg picture. For a 
two-level system, the Schrodinger equation is 



ih 



where ip has two components ipi ^^^d ijj2- The 
Hamiltonian can be written H{t) — {h'y /2){Bo{t) + 
J2j Bj{t)aj), where the constant 7 has been added for 
later convenience and where aj are the Pauli matri- 
ces, so that CTaCTb = 6ab + KahcC^c- Defining f{t) = 
exp(— 17 i3o(T)dr) and ip{t) — f{t)ilj'{t) turn the 
Schrodinger equation for ^ into a Schrodinger equation 
for ip' with the Hamiltonian H = {hy/2) BjCj. Thus, 
without loss of generality, we use the latter Hamiltonian. 

Following Feynman and coll. the relation between 
the density matrix p (with matrix elements pij = i^ii^*) 
and the magnetic moment M is = P12 + P21, My — 
i{pi2 — P21) and Mz = pu — P22- Thus, we obtain p — 
(1/2)(1 + J2jMj(Jj) and \M\ = 1. We see that £, = 
{Mx + iMy)/{\ + Mz) = '02/'0i- If we diagonalize p we 
recover the states t/'i and V'2 up to a phase that cannot 
be specified easily. Therefore, we shall work with the 
density matrix. The equation of motion for p is 



dp 

dt 



The commutation relations for the Pauli matrices turn 
this equation into the Bloch equation 

M = -7M X B. 

To turn the Bloch equation into the LLG equation, we 
just replace M by M' and B by B - (a/7)M'. If we 
denote by p' the density matrix corresponding to M', we 
find the equation of motion in the presence of damping 



P 



--[i7,p']+ia[p',p']. 



In other words, the damping term of the LLG equation is 
transformed into a nonlinear term ia\p' .p'\ in the equa- 
tion of motion of the density matrix. 
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If we replace B by B — (a/7)M' in the Schrodinger 
equation itself, we obtain the nonlinear Schrodinger equa- 
tion 



and the corresponding solution of the LLG equation is 



dV4 
dt 

dV4 
dt 



((1 + 2ia|^^p)i?3 - zaV'i(V'^)*i3+)^; 
-((1 + 2za|V;i')i?3 + za(^;)*Vii?-)V2 



Therefore, our method transforms the analytic solution 
of a Bloch equation into the analytic solution of a non- 
linear Schrodinger equation. We also reach the surpris- 
ing conclusion that the nonlinear terms in the nonlin- 
ear Schrodinger equation can describe the influence of 
damping. However, this damping does not create de- 
coherence: it transforms M into M', which is real and 
satisfies |M'| = 1. As a consequence, (p')^ = P' ^^'^ P' is 
the density matrix of a pure state. 

We investigate now the effect of damping on two fa- 
mous nonperturbative effects in two-level systems: self- 
induced transparency and dynamical localization. 

McCall and Hahn [l^ discovered a solution for the hy- 
perbolic secant pulse B = (a/ cosh(t/T), 0, 0). We recall 
that /^dt/cosh(t/T) = TTT. We consider more gen- 
erally a spin system submitted to a time-varying mag- 
netic field linearly polarized along Ox: B = (&(i),0,0), 
where h(t) is only required to be integrable. At time 
t = the spin has the spherical coordinates 6*0 and 0o- 
Let ]{t) = \l dr&(T) and a = (1 - eo)/(l + fo) with 
^0 ~ tan(0o/2)e*"*°. The solution of the Bloch equation 
is 



My{t) 

M,{t) 



M- 



psin(x(t) -I- u) 



-2M 



2M 



pcos{x{t) + u) 



with p = \a\, a = pe*" and x{t) — jfit). Self- induced 
transparency occurs when a;(oo) = 2mr because, after a 
long interaction with the external field, the system finds 
itself in its state at to. In the case of the McCall and 
Hahn pulse, we find f{t) — 2aTarctan( tanh((t— to)/2T)) , 
which tends to aTTr/2 for large t. Therefore, self-induced 
transparency occurs when jar — 4n for some integer n. 

To determine the effect of damping on this phe- 
nomenon, we calculate 



ae 



M 



^2ax{t) 



P 



M'Jt) = -2M 



(*^psin(x(t) -I- u) 



M'At) 



2M 



3"*(*)pcos(i(t) +m) 



Q2ax{t) _,_ p2 

with x{t) = 7/(i)/(l + a'^). Two effects can be observed. 
Firstly, the resonance condition is shifted from a::(oo) = 
2mr to a::(oo) = 2mr(l + a^); secondly, even at resonance 
the initial state is not fully recovered. For example, if 
the system is initially in the state ipi{0) — I, tp2{0) — 0, 
we have M — 1, p = 1 and w = and the final state is, 
at resonance. 



K(oo) 
M'^oo) 



= tanh(aa:(oo)), 
= 0, 

I 

cosh(aa;(oo)) 



Thus, it is not possible to recover the initial state because 
ax{oo) 7^ 0. Therefore, although the damping of the LLG 
equation amounts only to a trend towards an equilibrium 
state and not to a decoherence, it leads to a loss of self- 
induced transparency. However, if a is small, this loss is 
reasonably small and the reader might have the feeling 
that damping has only a minor effect on the dynamics. 
The next example shows that this is not the case. 

Quantum systems under the influence of a periodic or 
quasiperiodic external field can exhibit a freezing of the 
diffusion of its quantum state (l6j . called dynamical lo- 
calization. This effect is explained by the Floquet struc- 
ture of its quantum state. It was analyzed rigorously 
for a two-level system with B = (— /(t),0, e) by Barata 
and Cortez In the case where /(<) = acosuit, they 
showed that the crucial parameter is (g^), the time av- 
erage of q^{t), where q{t) — exp (ia7sin(ijji)/a;). Dy- 
namical localization occurs when (q^) = Jo(2a7/w) is 
zero, i.e. when x = 2aj/uj is a zero of the Bessel 
function Jq. In that case, dynamic localization is de- 
scribed as a modulation with a very slow secular fre- 
quency n = 2-f^e^T{x)/uJ^ + OCe-*), with T{x) = 
-J2m,n'^"ix)Jn-mix)Jmix)/i'mn), where J„ are Bessel 
functions and the sum is over all nonzero integers. 

In the presence of damping, 7 is replaced by 7/(1 — ia) 
and (2^) can no longer be zero because all zeros of Jq are 
real [T^l- Thus, contrary to the case of self- induced trans- 
parency, no resonance is strictly possible in the presence 
of damping. However, a second more dramatic effect oc- 
curs. If we take x = 2aj/uj a zero of the Bessel function 

X/(1 + a^), the Bessel 



Jo, and expa,nd (g^) around A 
expansion [I9| yields 

= Jo(A(I + za)) 



00 

fc=i 



(-l)'=Jfe(A)Jfe(iAa). 
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Therefore, 

= -zxJi(x)a-(x/2)V2(x)a' + 0(«'). 

The off- resonance secular frequency is ^ — ^liQ^}, 
so that 

n = ~ie^axJi{x)-^lo?{xl2fUx)+0{a^). 

The factor Ji (x) is not zero because the Bessel functions 
Jo and J I have no common zeros [l8| . In other words, 
the secular frequency acquires an imaginary part. Barata 
analyzed the influence of an imaginary part in detail [l^l 
and showed that the Floquet structure of the wavefunc- 
tion generates exponentially increasing terms that de- 
stroy the unitarity of the quantum evolution. We have 
shown that the norm of the wavefunction is conserved by 
the damped dynamics. Thus, the quantum evolution is 
still unitary and damping destroys the Floquet structure 
of the wavefunction. In this last example, damping has a 



dramatic effect on the dynamics of the system. Indeed, 
experiments confirm that dynamical localization is very 
sensitive to external perturbations (2lj |. 

It was shown in this paper that analytic solutions of 
the Bloch equations can be transformed into analytic so- 
lutions of the LLG equation. The influence of damping 
leads to interesting results for two-level systems, in par- 
ticular the transformation of the linear into the non linear 
Schrodinger equation. The effect of damping was shown 
to be quite different for self-induced transparency, where 
it is moderate, and for dynamical localization, where it 
destroys the structure of the dynamics. 
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